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ABSTRACT 


There are several issues that are encountered when designing a Free Electron Laser (FEL) 
oscillator in the terahertz (THz) regime. One such issue is that there are few inexpensive 
materials that are suitable for forming a semi-transparent out-coupling mirror. In an 
attempt to mitigate this particular issue, the focus of this thesis is to explore various out- 
coupling techniques for a notional FEL oscillator in the THz regime. The advantages and 
disadvantages of several out-coupling methods will be discussed, including hole out- 
coupling, and double mesh out-coupling. Simulation results will be presented, showing 
the effects of these various methods on FEL extraction and optical beam quality. 
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I. 


INTRODUCTION 


A Free Electron Laser (FEL) that uses electrons that are not bound to an atom to 
create coherent light. This type of laser has many attributes that distinguish it from the 
myriad of different lasers used worldwide. One characteristic is that FELs are tunable 
over a wide range of the electromagnetic spectrum and can be designed to create virtually 
any wavelength of light, from microwave to X-ray. This is because FELs do not depend 
on a conventional gain medium to produce the coherent photons; rather, the photons 
produced are a result of the electron interaction with a magnetic and optical field. The 
FEL can be operated at much higher energies than other conventional lasers since there is 
no gain medium to be damaged. 

This thesis will be exploring a notional FEL operating in the Terahertz (THz) 
regime of the electromagnetic spectrum. Operating in this regime does present some 
unique issues. The first issue presented is that of clipping, which will be explained in a 
later chapter. Briefly, clipping is an issue because of long wavelength associated with 
THz. Another issue is that many materials are not transparent to THz light. This makes 
it more difficult to extract usable light created in the FEL cavity. This issue is the 
motivation behind the study of various out-coupling techniques for an FEL. 

The methods of out-coupling to be explored will consist of hole out-coupling and 
mesh out-coupling. Hole out-coupling utilizes a highly reflective material with a hole in 
the center as the out-coupling mirror and mesh out-coupling uses a fine metallic mesh 
grid structure in place of the mirror. 

The notional FEL will be in the oscillator configuration. The components of this 
configuration will be described later in this chapter. 

A. FEL COMPONENTS 

A FEL in the oscillator configuration is comprised of several components. The 
major components as shown in Figure 1 include the injector, the accelerator, the 
undulator, the optical cavity, the energy recovery area, and the beam dump. 
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Figure 1. FEL in the oscillator configuration. Red dots are electrons, 
purple ellipses are light, and light blue ellipses 
are out coupled light. From [1] 


1. Injector 

The injector is where the “free” electrons are created from the cathode. The 
processes by which electrons are emitted from the cathode surface include thermionic 
emission and the photoelectric effect. Thermionic emission utilizes a heated cathode that 
emits electrons from its surface. The photoelectric effect uses light pulses that are 
incident on the photocathode. If the energy of the light is greater than the work function 
of the cathode, then electron pulses will be emitted [2], The emitted electron pulses are 
then accelerated using a static or oscillating electric field to around 5 MeV. 

2. Accelerator 

The electron pulse next enters the accelerator. Many accelerators are 
superconducting radio frequency linear accelerators (RF LINAC) that have niobium 
cavities that require liquid helium for cooling. In this stage, an electron pulse travels 
through each cavity in phase with the oscillating RF field inside the cavity. As a result, 
the electrons can be accelerated to a highly relativistic energy, 100 MeV for typical 
applications. The FEL can have electrons that have energies ranging from 5 MeV to 15 
GeV. 
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3. 


Undulator 


The undulator is composed of a series of permanent magnets that are arranged in 
alternating polarity and are often spaced equally. As the electrons pass through, the 
Lorentz Force causes them to “wiggle” as they travel. As the electrons wiggle, they 
create a light pulse that has a wavelength much shorter than the undulator period because 
of relativistic Doppler shifts. Only a small percentage of the energy of the electron pulse 
is converted into light. 

4. Optical Cavity 

The optical cavity is the region of the FEL where coherent light is amplified. The 
components of this region are only used in the oscillator configuration. The oscillator 
configuration uses a short undulator and two mirrors. One of the mirrors is highly 
reflective and reflects the created light pulse back through the undulator for stimulated 
emission and coherence evolution. The out-coupling mirror is typically partially 
transparent and allows for a small percentage of the light to pass through. 

5. Energy Recovery (The Beam Dump) 

Since only a small percentage of the energy of the electrons is converted to light, 
some FELs direct the electrons to pass through the RF accelerator again with the RF field 
~ 180-degrees out of phase. This causes the electrons to lose energy to the RF field, and 
prevents high energy electrons from entering the electron beam dump. This reduces the 
amount of shielding needed for any applications. This type of configuration is call an 
Energy Recovery LINAC (ERL), since it recovers most of the unused energy in the 
electron pulse, thus improving the overall efficiency. 
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B. FEL THEORY 

1. FEL “Pendulum” Equation and Electron Dynamics 


The dynamics of electron motion as it passes through the undulator are governed 
by the following equations: 


d{yp) 

dt 


e 

me 


(.E + fixB ) 


dy 

dt 



me 


( 1 . 1 ) 


( 1 . 2 ) 


r = 



(1.3) 


where E and B are the electric and magnetic fields, respectively, e = \e\ is the charge 
magnitude of an electron, m is the mass of an electron, /3 = v/c is the dimensionless 

electron velocity, and y is the relativistic Lorentz factor [1], Equation 1.1 is the Lorentz 
force equation and governs the motion of the electron as it encounters FEL fields. 
Equation 1.2 is the rate of energy change due to the work done by the electric field. This 
equation calculates the energy given to the optical field by the electron. 

As the electrons travel through the undulator, they interact with the electric and 
magnetic fields of the light, E = £'(cos^,-sin^,0), B = E (sin y/, cos ^,0) where E is the 
field amplitude in cgs units, y/= kz-cdt + tf), with (f) the optical phase, k = 2njA the 
wave number, A the wavelength, and co the optical frequency. The electrons also 
interact with the magnetic field of the undulator, B = 5(cos(£ 0 z),sin(& 0 z),o), where B is 
the undulator field strength, and A 0 = 2 n/k^ is the undulator period [1]. The electrons at 
this point are highly relativistic, so fi z ~ 1, so when the fields are substituted into 
Equation 1.1, the transverse motion in the undulator is approximately determined only by, 
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(1.4) 


J3 ± «-(cos(& 0 z),sin(k 0 z),o) 


r 

9 

where K - eB iws A () j^27imc 2 ) is the dimensionless undulator parameter, and B rms is the 

rms value of the undulator field. Typically K « 1 for most FELs [3]. If Equation 1.4 is 
substituted into Equation 1.2, the electron’s energy evolution yields 

eKE 

cos + 

(1.5) 


y = — = —cos^-/? sin^~| = ^^cos(<^+^) 
dt me L y J ymc v ’ 


where =(k + k^z-cotis the “electron phase.” The electron phase follows the 
microscopic position of the electron as it evolves. It is important to note that the energy 
exchange between the electron and the optical field increases when cos(<^+ ^) > 0 and 

decreases when cos(^ + (f) < 0. 


The “electron phase velocity” is defined as v = —£=L[[k + k 0 )/3,-k^ = £, 

where (...) indicates a derivative with respect to dimensionless time t = ct/L with L the 

length of the undulator so that T varies from 0 to 1 [1], If v= 0, the energy exchange is 
at an optimum level and the FEL is said to be at “resonance.” In this condition, the 
electron phase velocity is manipulated to show that the resonant laser wavelength is 


t 4>(i-A) ., 4>(i-* 2 ) 
A ~ 2 7o 2 


( 1 . 6 ) 


where yo is the initial electron beam Lorentz factor and y 0 » 1. 


In a beam of electrons, each electron will have a unique initial phase 
do )=(«=( k+ Kh as they enter the undulator. In a typical FEL, ~10 6 electrons are 

randomly spread in each section of the electron beam ~/lin length. This means that 
roughly half of the electrons will lose energy to the laser field and half will take energy, 
as described in Equation 1.5. A change in the electrons energy will also create a change in 

the position, ^and the electron phase velocity, v. The change in the phase velocity is 
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i>= 4;ry/y 0 while near resonance in the relativistic limit [1], If this is combined with 
Equation 1.5, then the FEL “pendulum equation” is derived 


v = 4' = Hcos(4'+^) 5 ( L7 ) 

where |a| = 4 ftNeKLE / y^mc 2 is the dimensionless laser field amplitude[l]. From this, it 

is shown that changes in the design of the FEL will affect the dimensionless amplitude of 
the laser field. If the electron beam has a larger Lorenz factor, the strength of the laser 
field will need to increase in order to create the same dimensionless amplitude. 

2. FEL Wave Equation 

Having solved for the microscopic motion of the electrons in the undulator, the 
evolution of the laser field in the undulator must be derived. Consider the full vector 
wave equation with the Coloumb gauge [3] 


V- 


1 d 2 ^ 
c 2 dt 2 . 




( 1 . 8 ) 


where J± = —/? £(3) 

[x-rt (t)j is the transverse current density, r, (?) is the electron 

i 

position at time t, and A^x,t ) is the optical vector potential. The electric and magnetic 
fields of the laser can be derived using the vector potential, 


E = -— _ 

c dt B = V x A 


(1.9) 


For the case of laser light, the beam is partially or completely coherent and 
A (x,t j has a slowly-varying envelope E(x,t ) in the direction of propagation along the z 

axis [3]. If the electric field is assumed to be complex and slow-varying over an optical 
wavelength, the optical vector potential is assumed to have the form 


A 


^x,t) = ^^-se i[kz ~ at) 


(1-10) 
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where E = \E\e" /> is the complex laser electric field, and £is the polarization vector of the 




and the phase 


laser field. If the electric field amplitude 
be slow-varying along the z axis [3], Equation 1.8 can be rewritten as 

i(kz-6)t) 


x,t) 


are assumed to 




—*2 ^ . T 1 

r d 15) 

Vi + 2ik\ 



v dz c dt) 


An - 
E*—-Ji 
c 


an) 


Both sides of the equation are multiplied by ke ,{t <0 ^£ and a new coordinate, 
u = z-ct, is introduced in order to invoke the method of characteristics!!]. With that, 
Equation I! 1 becomes the paraxial wave equation with a source current 


—-2 ^ i 

n dY 

Vi + 2 ik\ 



vc dt) 


E = -—J 1 .£* e - i{kz - a,l) 


( 1 - 12 ) 


Using the relationship between the transverse motion of an electron and the 
transverse current density defined previously, Equation I! 2 becomes 


—2 2 ik d 

v±+—— 

c dt 


E = - AmeKkp{x,t ) < e ^ jy 


> (P) 


(1-13) 


where pyx,tj is the local electron density in a small volume element dV located at yx,tj, 
and <.. .> represents an average over sample electrons in that volume element. 

Dimensionless time, t, is reintroduced and the wave equation is multiplied by the 
factor -AnNeKL 2 /ylmc 2 k so that it becomes 


iL -2 d 

-V_|_ H- 

2k dr 




(1.14) 


where a = \a\e" f ' is the complex laser field, \ci\ = AjiNeKLEjypnc 1 is again the 

dimensionless laser field amplitude, j = $7r 2 Ne 2 K 2 L 2 pjy 3 0 mc 2 is the “dimensionless FEL 
density,” and it is assumed that y ~y 0 for all the electrons in the beam[3]. 
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The wave equation is made completely dimensionless by defining the 
dimensionless transverse coordinates x ~ x {k/2L) , anc [ >' The 

equation is now written as 


wave 


l =2 


--V1 + — 
4 dr 




(1.15) 


where the coordinates (x, y, z, r) are now all dimensionless. When the laser beam and the 
electron beam overlap and diffraction is small over the undulator length, the FEL wave 
equation can be written in its simplest form 


a = -j <e ^ > 


(L 16 ) 
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II. THZ ISSUES 


A. MATERIAL CONSIDERATIONS 

For the purpose of this thesis THz radiation is defined to be in the range of 100 
microns to 1mm. There are several issues that arise when out-coupling THz radiation in 
a FEL oscillator in this range. Typically, FELs use a partially transparent mirror to out- 
couple the light from inside the optical cavity. These mirrors consist of a substrate that is 
then layered with a coating that is reflective to the wavelength of light being used. As 
such, there are several material properties of the mirror substrate that must be accounted 
for in order to produce usable THz light for the FEL. 

1. Transmission of Radiation 

In order for a FEL to operate, a fraction of the light produced in the optical cavity 
must be transmitted through the out-coupling mirror of the resonator. THz radiation is 
absorbed by many materials, which often leads to the need for expensive materials to be 
used. 

The fraction of incident radiation absorbed can be calculated by: 

A = \-{\-R)e~ ax 

where A is the fraction of radiation absorbed, R is the reflection coefficient, a is the 
absorption coefficient, and x is the thickness of the material [4], This shows that with an 
increase of the thickness of the material or with an increase in the absorption coefficient, 
the fraction of light absorbed increases. This absorbed energy manifest itself as heat and 
can result in other effects on the mirror. 


2. Thermal Expansion 

When a majority of materials are heated by a source, there is a tendency for the 
materials to expand. This thermal expansion is found experimentally and calculated by: 


a T 


L 


dL\ 
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where aj and f> are the coefficients of linear and volumetric expansion, respectively, L 
and V are the materials length and volume, respectively, and T is temperature [5]. From 
these two equations, it can be seen that in order to have a material that exhibits minimal 
change in dimension, the material in question must have a small coefficient of expansion. 
A minimal change in dimension is critical in order to prevent optical distortions, which 
can lead to the laser either not operating at the optimal levels or not operating at all. 


3. Thermal Conductivity 

Thermal conductivity describes the ability of materials to conduct and transfer 
heat. Thermal conductivity is found experimentally and calculated by 

k= L fdQ' 

dt j 

where k is the thermal conductivity in watts per meter-Kelvin, L is the lateral distance 

between the temperature difference AT, and is the rate of heat flow into the 

material. [6] The thermal conductivity is shown to be proportional to the heat flow. As 
such, a material with a large k value will be able to remove heat better than one with 
small values. 

Because the out-coupling mirror absorbs large amounts of energy from the FEL, it 
is important that the material that comprises the lens be able to transfer the heat from the 
focal point of the beam. If this is not achieved the beam will distort and/or potentially 
destroy the mirror. 


4. Insufficiencies 

There have been some studies and experiments conducted to identify the optimal 
material for use with THz radiation. Table 1 provides a summary of selected materials 
with experimental data for the previously discussed properties. The color scheme of the 
table is used to highlight material property values in a range that would be suitable for 
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use in a FEL operating in the THz range. As shown, only one of the five materials, 
Chemical Vapor Deposition (CVD) Diamond, appears to be suitable for use as a mirror. 
However, even CVD Diamond has its drawbacks. Because CVD Diamond is very 
expensive, this would lead to a very expensive mirror regardless of the dimension. 


Materials 


Thermal 
Conductivity 
@20 °C (W/M °C) 


Thermal 
Expansion 
Coefficient 
(xlO-6 / °C) 


Transmission 

Range 

(pm) 


HRFZ Si 

1.5 

2.6 

50-300 

Quartz 

1.3 

0.59 


Sapphire 

42 

5 


Zinc Selenide 

18 

7.57 

10-20 

CVD Diamond 

3300 

1.2 

10 - 300 


Table 1. Summary of selected materials and respective properties. From [7] 


B. CLIPPING 

FEL Beam Clipping is the process by which an optical beam is negated when it 
comes in contact with an absorbing material. Because of this interaction, a portion of the 
beam’s energy is lost which results in a decrease in power. This phenomenon is most 
prevalent in the undulator section of the FEL because the beam line is constricted. 
Figure 2 illustrates how a clipping interaction in the undulator would occur. 



Figure 2. Clipping in the undulator as a function of diffraction. 


11 




























A contributor to the degree of severity of the clipping that occurs is beam 
diffraction. Beam diffraction can be described by 


w(z) = w 0 



(II. 1) 


where w(z) is the radius of the beam at distance z, wo is the optical mode radius and z R is 
the Rayleigh length [1], From Equation II. 1, as the distance from z=0 increases the radius 
of the beam also increases. This effect is further exaggerated because the Rayleigh 
length, defined by 


n nw; 


Z R=- 


0 


/l 


(II.2) 


where A is the optical wavelength, and n ~1 is the index of refraction decreases. As the 

wavelength increases the Rayleigh length increases which in turn increases the beam 
radius at distance z. This causes the beam to expand rapidly and as such increases the 
amount of the beam that potentially will be clipped. 

The potential for beam clipping is very prevalent with THz radiation and must be 
taken into account in accordance with previous studies [7] .As such; the simulations that 
will be conduct will use a constant clipping factor of 4% per pass through the undulator. 
A separate study of how clipping affects the extraction of energy of the FEL operating in 
the THz regime will also be conducted in a later chapter. 
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III. GAUSSIAN MODES 


A. GAUSSIAN BEAMS 

There are several ways to model the light wave produced by an FEL oscillator. 
The most convenient approach is to use Gaussian modes. A Gaussian beam in the 
fundamental transverse mode describes a light wave where the power and intensity of the 
light are at a maximum at the center of the beam and decrease as the distance from the 
center increases. This approximation is a solution to the paraxial wave equation. 

1. Short Derivation of Paraxial Wave Equation 

Consider the wave equation found by solving Maxwell’s equations in free space: 

v2 £=« 02 % 2 ; (mi) 

where the spatial derivatives are taken to be in rectangular coordinates. If it is assumed 

that E is polarized in one direction, so that the scalar form of the wave equation can be 
invoked 

= (1,1.2) 

where l/ c 2 = ju () s () . It can be seen that E has a spatial and temporal dependence. As 

such, the electric field can be defined as the product of both spatial and temporal factor so 
that 

E(x, y, z, t) = E{r , t) = E(r)e iat ; (III.3) 

where co = (2/re) j X = kc . When substituted into III.2, the wave equation becomes 

[v 2 + k 2 )E = 0; (III.4) 

where k 2 = of j c 1 ,and the explicit time dependence has been removed from the wave 
equation [8]. 
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To arrive at the paraxial wave equation from Equation III.4, it is assumed that the 
transverse beam profile changes slowly on the scale of an optical wavelength. As such, 
we assume the form 

E(x,y,z) = u(x,y,z)e~ ,h -, (III.5) 


where u describes the transverse profile of the traveling beam [8]. In other words, u 

changes slowly in z relative to the factor, e ~ ,kz . When Equation III. 5 is substituted into 
III.4 the wave equation becomes 


a 2 ~ d 2 ~ d 2 ~ d ~ 

— -u-\ - -u-\ - -u-2ik—u = 0. 

dx 2 dy 2 dz 2 dz 


(III.6) 


In III.6 the second derivative of the z component of u can be assumed to be negligible by 
invoking the paraxial approximation, 3 2 «/<9z 2 | <§: \2kdu/dz\ [8], From this, Equation III.6 
can be rewritten as the paraxial wave equation: 


f 


d 


V -2 ik—, 
1 dz) 


u = 0; 


(in. 7) 


where V 2 = d 2 /dx 2 + d 2 /dy 2 . This approximation remains valid for traveling beams that 

propagate at an angle less than 30 degrees with the z-axis. Due to this, the paraxial 
approximation holds true for most FELs in the oscillator configuration. 

B. HERMITE-GAUSSIAN BEAMS 

A solution to the paraxial wave equation must be found in order to properly model 
a propagating beam. One such solution is the Hermite-Gaussian solution, 


E lm (x,y,z) = E 0 -^H l 
w(z) 


( \ 
w{z) J 


H„ 


f \ 

V2 y 


v 


w(z) 


) 


xexp 


-ik - + ~^ - X —ikz + i(l + m + l)^(z) 


27? (z) w 2 (z) 


(III.8) 
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This solution is very intricate yet it is still in the form of E\x,y,z ) = u[x,y,z^e~ lkz 
[8], The waist size term, w(z), was initially defined by Equation II.l. R(z) is the radius of 
curvature of the propagating beam and it is given by R[z} = z[\ + [z 2 R / z 2 ^. From this 

definition, it can be seen that the radius of curvature changes depending on whether the 
beam is propagating in the positive or negative direction. rj(z) represents a phase shift 

defined by ^^zj = arctan^z/z fl j. The / and m labels are integer indices of the Hermite 

polynomials, Hi and H m . These indices ultimately affect the transverse intensity of the 
beam and are referred to as modes. The 1=0, m=0 mode is called the “fundamental” 
mode because the subsequent Hermite polynomials are both more structured. Figure 3 
displays several of the lower order spatial modes that the beam can assume. 



Figure 3. Transverse Electromagnetic Modes. From [9] 

For the purpose of the FEL, the fundamental mode is often preferred because the 
beam can be later focused into a tighter spot for a given set of focusing elements. For the 
purposes of directed energy, this means more power can be delivered to a smaller spot on 
a distant target. 
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IV. HOLE OUT-COUPLING 


With hole out-coupling, the light created in the FEL optical cavity is transmitted 
through a hole in a highly-reflecting mirror instead of being transmitted through a 
partially transparent mirror. This technique has been shown to work for other types of 
lasers operating at various wavelengths [10]. 

A. DESIGN PARAMETERS 

For the hole out-coupling simulations, Table 2 lists the design parameters that 
were held constant: 


E b ~ 5 MeV 

f ~ 50 MHz 

N= 16 

Xq ~ 2 cm 

X ~ 100 pm 

K ~ 0.35 


Table 2. FEL Parameters 


where Eb is the beam kinetic energy, f is the pulse repetition frequency, N is the number 
of undulator periods, Xq is the undulator period, X is the optical wavelength, and K is the 
dimensionless undulator parameter.[3] 

B. SIMULATION RESULTS 

Simulations were conducted using the NPS 3D oscillator FEL code and were 
categorized into two groups based on the electron bunch charge, lOOpC and the 200pC. 

1. lOOpC Bunch Charge Simulations 

Simulations were conducted for the lOOpC bunch charge varying the radius of the 
hole, the quality factor Q, and the clipping gap. 

a. Varying Hole Radius with a high 0 value 

In the first set of simulations conducted for the lOOpC bunch charge, the 
hole radius was varied over a range of values. The quality factor, Q, is defined to be 
inversely proportional to the amount of loss in a cavity without a hole. For this set of 
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simulations, Q was set at a very large value (~10 10 ) in order to approximate a highly- 
reflecting mirror surface. This means that essentially no light is absorbed by the mirrors, 
thus the only way for the light to escape is through the hole. The hole radius range was 
from 0.05 cm to 0.40 cm. The simulation results are plotted in Figure 4. The shape of 
the plot shows that there is a peak of -2.8% extraction at a hole radius of 0.125 cm. This 
is due to how the electrons are evolving to produce the light, as will be seen in 
subsequent figures. 


Extraction for a Highly Reflective Surface with a Hole 



Figure 4. Simulations results for lOOpC bunch charge with Q= 10 10 

Figure 5 (a) displays the optical field evolution in the undulator for the 
hole radius of 0.05 cm. The blue/cyan area represents the optical field, the red dots are 
the electrons, and the yellow curves represent the mode of the beam at the ends of the 
undulator. In this case the peak dimensionless optical amplitude is high, |a|«52. 

indicating FEL saturation in strong optical fields.. With a small hole only a small 
fraction of the light produced in the optical cavity is out-coupled. The optical field grows 
to saturation, while the electrons also evolve described by the “pendulum” equation and 
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the FEL wave equation. In Figure 5 (b), the electrons appear to overbunch as they travel 
farther in their phase space orbits, absorbing energy from the optical field and hence 
reducing the extraction. 



Figure 5. Simulation output for hole radius of 0.05 cm. (a) The optical field evolution in 
the undulator. (b) The electron evolution phase space plot at the end of the undulator. 


The optimal hole radius of 0.125 cm gives the largest amount of optical 
power out of the FEL mirror hole. The peak optical field is smaller (|a| = 27) than that of 

the smaller hole radius, since a greater fraction of the light is out-coupled per pass. Figure 
6 (b) shows the electrons do not overbunch as severely, and hence more energy is 
converted to light and out-coupled through the mirror. 



0 t 1 


(a) (b) 

Figure 6. Simulation output for hole radius of 0.125 cm. 

If the hole radius becomes too large, too much of the light escapes from 
the cavity which causes a drop in the dimensionless field amplitude (|a| = 2.4) as seen in 
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Figure 7 (a). Due to this, the electrons do not evolve much in phase space, as seen in 
Figure 7 (b), and very little energy is extracted from them. 



(a) (b) 

Figure 7. Simulation output for hole radius of 0.25 cm. 


The use of hole out-coupling did have a negative effect on the transverse 
mode shape. As seen in Figure 8, as the size of the hole increased, small perturbations in 
the mode shape at the out-coupling mirror started to develop. The beam is losing energy 
to the formation of higher-order modes, which is not a favorable condition for the FEL 
interaction and applications. 


a (x, y, 5 . 5) 6 . l|a (x, y, 5 . 5) 3 .6] a (x, y , 5 . 5)0 . 3 

^ ^ 


(a) (b) (c) 

Figure 8. Transverse mode shape at the out-coupling mirror for a hole radius of:(a) 0.05 cm, 

(b) 0.125 cm, and (c) 0.25 cm. 

b. Varying Hole Radius with a Q value of 50 

For this set of simulations, the quality factor was Q= 50, representing a 
loss of 1% of the energy due to absorption at each mirror surface. The hole radius was 
again varied from 0.05 cm to 0.40 cm. The range was selected to provide a direct 
comparison between the highly-reflecting mirror surface and the more realistic mirror 
surface. 
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Figure 9 plots the overall results of the simulations showing a trend 
towards smaller hole radii. A smaller hole radius is expected to be necessary, since now 
that some of the energy is absorbed by the mirror surfaces, the out-coupling hole should 
be smaller to compensate for the absorption. Due to limitations on the numerical grid, the 
optimal hole radius was unable to be found, which appears to be less than 0.05 cm for 
this case. 


Extraction for a Partially Absorbant Surface with a Hole 



Figure 9. Simulations results for lOOpC bunch charge with Q= 50 

c. Varying the clipping gap for the theoretical simulations 

Since the longer wavelengths expand more rapidly than shorter 
wavelengths, absorbers were simulated at the ends of the undulator to clip a portion of the 
field that might scrape in the ends of the undulator. The gap between the absorbers was 
varied from 0.3 cm to 3.0 cm while the hole size was kept constant at 0.125 cm. For the 
results in the previous section, the absorbing gap was kept constant at 1.4 cm. 

Figure 10 shows the results of the simulations, and finds that there is an 
optimal clipping gap of 1.5 cm for this set of parameters. As the gap size increases past 
the optimal there is little change in the extraction levels. This intuitively makes sense 
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because if the size of the gap keeps increasing, eventually there is no clipping. Past that 
threshold, there should be no change to the extraction since no energy is being lost to the 
absorbers. 


Gap Varitions for Bunch Charge of 100 pC 



Figure 10. Simulation results for Gap variations. 


Figure 11 shows the simulation output for the optimal gap of 1.5 cm. It 
can be seen in that the mode shape is starting to become slightly less Gaussian in (a) and 
(c). This means that slightly more energy is going into higher-order modes, and shows 
that clipping can be used to force the optical mode into the fundamental. 



Figure 11. Simulation output for gap of 1.5 cm. 
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2. 200 pC Bunch Charge Simulations 

Simulations were also conducted for the 200 pC bunch charge, again varying the 
same parameters: the quality factor Q, the radius of the hole, and the clipping gap. 

a. Varying Hole Radius with a high Q value 

In the same way as with the 100 pC simulations, the set of simulations 
were conducted over a range of hole radii with Q ~10 10 . The same hole radius range of 
0.05 cm to 0.40 cm was utilized. 

The results of the simulations are shown in Figure 12. The results show 
an optimal value for the hole radius of 0.15 cm with a corresponding extraction of 2.9%. 
The optimal hole size is slightly larger than the 100 pC case since there is more gain per 
pass with the 200 pC bunch, which allows more losses in steady-state. 

Extraction for a Highly Reflective Surface with a Hole 



Figure 12. Simulations results for 200 pC bunch charge with Q= 10 10 

Figure 13 displays the output results for the optimal hole radius of 0.15 cm 
showing a well-bunched distribution in phase space in Figure 13 (b) and slight 
perturbations in the mode shape due to the hole in Figure 13 (c). 
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b. Varying Hole Radius with a Q value of 50. 

For this set of simulations, the quality factor was Q = 50. As previously 
mentioned, this represented a loss due to absorption at the mirror surfaces. The hole 
radius range again was from 0.05 cm to 0.40 cm. 

Figure 14 illustrates the overall results of this set of simulations directly 
compared to the 100 pC case. It shows that there was an increase in the optimal hole 
radius compare to the smaller value as seen in the previous case. For this set of 
simulations, the optimal hole radius is 0.10 cm and produced an extraction of 2.9%. The 
optimal hole radius is larger than the 100 pC case, since the higher gain per pass from the 
200 pC bunch charge can tolerate more round-trip losses due to absorption and out- 
coupling. 

Extraction for a Partially Absorbant Surface with a Hole 



Figure 14. Simulations results for 200 pC bunch charge with Q =50 
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Figure 15 shows the a comparison between the optimal hole radius of 0.05 
cm for 100 pC case and the 200 pC optimal hole radius of 0.15 cm both with a g~50. 
Although the dimensionless amplitudes are not very different, the mode shape on the out- 
coupling mirror shows that this set of parameters will include higher-order modes as the 
hole size increases. 




Figure 15. Simulation output for hole radius of (a) 0.05 cm with 100 pC and (b) 0.15 cm with 

200 pC. 


c. Varying the clipping gap for the theoretical simulations 

Again, simulations were conducted varying the size of the absorbing gap 
over the range of 0.3 cm to 3.0cm. 

Figure 16 shows the results of the simulations. The optimal gap in this 
case is 1.2 cm with an extraction of 3.0%. After the peak extraction was reached, the plot 
leveled off much like it did for the lOOpC set. 


25 












Gap Varitions for Bunch Charge of 200pC 



Figure 16. Simulation results for Gap variations. 
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V. MESH OUT-COUPLING 


Mesh out-coupling is another out-coupling technique that does not use a 
conventional partially transparent mirror. Instead, the light in this technique will out- 
couple through a metallic mesh. There are two types of meshes that can be considered 
for this application. The first type is the typical metallic mesh that uses thin strips of 
metal with no substrate and is referred to as an inductive mesh [11]. The second type 
uses a thin metallic pattern deposited on a substrate and is known as a capacitive mesh 
[11], Both can have either circular or square symmetry and are illustrated in Figure 17. 
The mesh with square symmetry can be parameterized by three variables: the mesh 
period, g, the grid spacing, a, and the mesh thickness, t. 



Figure 17. Illustrations of freestanding (a) inductive metallic mesh and 
(b) capacitive mesh with square symmetry. 

From [11] t is the thickness of the mesh, g is the mesh period, 
and a is the grid spacing. 

Mesh out-coupling has been proven to work in the far infrared region of the 
spectrum [12]. While there are several ways to implement this technique [13], two will be 
discussed in this thesis: single mesh out-coupling and double mesh out-coupling. Both 
techniques will use an inductive mesh as the out-coupler in the thin mesh regime where 
t «a , with the assumption that there are no losses due to absorption in the mesh. 

A. OPTICAL PROPERTIES OF A SINGLE MESH 

The optical properties of a single mesh can be derived and analyzed to show that a 
single mesh can be utilized as a reflective surface. Previous study of this interaction has 
provided results that were used in creating parameters for this thesis [12]. The analysis 
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from this study showed that a single mesh of square symmetry has optical properties that 
are independent of polarization because the mesh can be considered as the superposition 
of two crossed strip gratings if the ratio a/g< 0.2is satisfied [12]. The same paper 
argued that if g/ A <1, i.e. when the wavelength is longer than the mesh period, 
diffraction from the mesh edges is negligible [12]. 


Many of the coefficients needed to describe the optical properties of the mesh can 
be calculated based on an equivalent electric circuit that yields: 

r = . 7/2 


1 + 7/2 
t= i+r 


(V.l) 
(V.2) 

where T is the complex amplitude reflection coefficient of the mesh, r t is the complex 
transmission coefficient of the mesh, Y - -i{2/Z^Aj g'j is the admittance of the 
equivalent circuit, and Z 0 = 21n(csc(;ra/g)is the impedance [12]. This assumes that 
there are no losses due to absorption in the mesh material. 


The theoretical power reflectance is the magnitude of the complex amplitude 
reflection coefficient of the mesh squared, R=\T\ . This can be used to show how the 
single mesh reacts to light that is normally incident onto it. In Figure 18, the power 
reflectance was plotted as a function of g/A and parameterized over 0.04 <a/g< 0.20 . 
The figure shows that when X is much larger than g, the single mesh tends to reflect most 
of the power. From this, it can be seen that the single mesh can be used in place of a 
mirror in the optical cavity of a long wavelength FEL, with the mesh parameters chosen 
to produce a desired quality factor for the cavity for a fixed wavelength. 
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Power Reflectance of a Single Mesh Out-Coupler 



Figure 18. Power reflectance of a single metallic mesh vs. g/A where each curve is 

parameterized by a/g. After [12] 


B. DOUBLE MESH OUT-COUPLING 

While a single mesh can be used as an out-coupler, a more flexible design uses a 
double mesh. Instead of using one freestanding metallic mesh, two freestanding metallic 
meshes are used in parallel and separated by some distance typically on the order of an 
optical wavelength, A. Then, the total transmission and reflection from both meshes 
depends not only on the properties of each individual mesh, but also by the separation 
between them, much like a Fabry-Perot Interferometer (which will be discussed in the 
following section of this chapter). Whereas the optical properties of the single mesh were 
strongly dependent on A, the extra degree of freedom afforded by a double mesh allows 
one to maintain the constant set of optical properties across a range of wavelengths by 
adjusting the separation between the meshes [14]. 
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1. Fabry-Perot interferometer 

A Fabry-Perot interferometer (FPI) consists of two partially transparent, parallel 
mirrors separated by a medium or vacuum. The FPI is also used in a variety of 
applications including but not limited to optical resonators and spectrometers. 

Figure 19 illustrates how a FPI transmits and reflects light, where s is the distance 
separating the two parallel reflecting surfaces, 6 is the internal angle of incidence, and E 0 
is the initial amplitude of the incident wave. As the wave travels through the FPI, each 
reflection creates a partial wave that is either transmitted or reflected. The subsequent 
internal reflections create, in principle, an infinite number of partial waves [15], 



Figure 19. Interaction between an incident light beam and a FPI. 


Each partial wave acquires a phase shift per pass through the interferometer 
according to 




Anns cos 6 

A 


(V.3) 


where n is the index of refraction of the medium between the two surfaces, and A is the 
wavelength of the incident wave, Eq,i n vacuum. The partial waves that are reflected (Eo, 
E 4 , Ee ,...) are then calculated as 


: rE {) , E 4 = t 2 rE 0 e ,s 


E 6 = t 2 r%e 2iS 
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where r and t are the amplitude reflection coefficient and the amplitude transmission 
coefficient respectively and are assumed to be equal for all reflections and transmissions. 
Similarly, the partial waves that are transmitted, {E\, £3, E 5 , ...), are also calculated as 

Ey = t 2 E 0 , E 3 = t 2 r 2 E 0 e iS , E 5 = t 2 r 4 E 0 e 2i \ 

The sum of these waves constitute the amplitude of the transmitted wave and the 
reflected wave, 



E t +E 3 +E 5 +...+E lcr _ x 

(V.4) 


E r =E 2 +E 4 +E 6 +...+E 2a 

(V.5) 

where E T and E R are the total transmitted amplitude and the total reflected amplitude, 

respectively, and a approaches infinity. When the values of the partial waves are 

substituted into Equation V.4 and Equation V.5, they become 


E T = t 2 (l + r 2 e iS + r 4 e 2iS + ..)E 0 , 

(V.6) 


E r = [r -+ t 2 re iS ( 1 + r 2 e iS + r 4 e 2iS + ...)]^ 0 . 

(V.7) 


It can now be seen that the total complex amplitudes both contain an infinite 
geometric series [15]. As such, the amplitudes can be rewritten as 


T 

Et " \^ E ° 


(V.8) 


(l 


1 -Re 1 


id 


F 

J-'n 5 


(V.9) 


where R = r 2 is the fraction of intensity reflected at each interface and T = t 2 is the 
fraction of intensity transmitted at each interface. 


If the intensity of initial wave is defined as I 0 - E* 0 E 0 , the ratio of transmitted 
intensity to the initial intensity can be described as 

I_t KEt l 1 "* 2 ) 

7 0 E* 0 E 0 (l-i^ + dtfsin 2 ^/!)' 
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(V.10) 



Similarly, the ratio of reflected intensity to the initial intensity can be described as 


T F* F 
1 r _ FrFr _ 


4Rsin 2 (8/2) 


(V.ll) 



(l-R 2 ) + 4Rsm 2 (S/2) 


An assumption of Equations V.8 through V.ll is r 2 +t 2 =1, i.e. there is no absorption 
occurring in the mesh material [15]. 

Figure 20 illustrates how the reflected energy depends on the phase shift for a 
given R. It can be seen that, even at large values of R, there are phase shifts where a large 
amount of the intensity is not reflected. From equation V.l 1, when 8 = 2m7i, where in is 
an integer, the reflected intensity is equal to zero. As such, these phase shifts allow for 
complete transmission, assuming that there is no absorption. 


Ratio of Reflected Energy to Incident Energy 



10 


5 


Figure 20. Reflectance of a single metallic mesh. After [15] 


Figure 21 displays the transmission ratio as a function of the phase shift. As seen, 
transmission maximums occur for the same phase value as reflection minimums. 
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Ratio of Transmitted Energy to Incident Energy 



Figure 21. Transmittance of a single metallic mesh. After [15] 

When using two meshes to form an FPI, the only change to the FPI equations is 
the reflectance, R, is equated to the square of T from Equation V.l. 

2. Spherical Double Meshes 

The NPS FEL 3D code assumes that the optical cavity is centered about the 
undulator and each mirror is assigned a radius of curvature based on the Rayleigh length 
of the cavity. This requirement motivated the consideration of spherical double meshes 
to ensure the feasibility of its use in the cavity [16]. 

Suppose a light beam has a normal incidence, 6 = 0, on the interferometer shown 
in Figure 19. In that case, the path length the light would travel through the 
interferometer is 2s and is the same for any beam normally incident along the 
interferometer. Therefore, all normally incident beams exit the interferometer with the 
same phase change. This ensures that the transmitted and reflected beams have the same 
phase. If a set of spherical surfaces are considered for the interferometer, it must be 
proven that the path length traveled is approximately uniform for any normally incident 
beam at any point on the surface of the interferometer in order to preserve the uniform 
phase change of the reflected and transmitted beams. 
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Consider two identical spherical meshes that are separated by a distance s, as in 
Figure 22. 



Figure 22. Beam interaction with spherical double mesh. Illustration is not to scale as z ~ 

meters, x ~ centimeters, and s ~ X. 

It must be proven that there is no appreciable change in the phase of the beam as 
it travels through both meshes as a function of the divergence angle, (p\ namely, 

Ac) = S(cp) - 8 (0)« 2n (V.12) 

where AS is the change in the phase. From Equation V.3, it is shown that, 
S ( ( p ) = (4 ns ') j A = (4tt,s)/(/1 cos (p) and S (0) = 4 tts/ X . When substituted into Equation 
V.12, the result is 

AS = [s/cos^-s]. 

Figure 22 shows that (p = tan 1 (x/z)«tan‘(10 2 ) «10 2 rad. From this, the small angle 

approximation can be invoked, cos cp « cos(10“ 2 )«1. Figure 22 also shows that s ~ X. 
The result is that AS «: 2n. 

3. Equating Quality factor, Q, to Mesh Spacing 

In order to simulate the out-coupling interaction of the double mesh out-coupler, 
the metal of the mesh was assumed to have negligible absorption for THz radiation. With 
this assumption, the only losses are transmission through the mesh. As such, the quality 
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factor, Q, that is defined to be inversely proportional to the loss of the cavity, is also 
inversely proportional to the energy transmitted through the double mesh, Q^T~ l . If 
the double mesh out-coupler is analyzed as an equivalent circuit with normalized 
admittances of each mesh being Y\ and Y 2 [12], the power transmittance of the double 
mesh out-coupler can also be calculated in the following fashion: 

rp I |2 

* -| r zw| ’ 

where z DM = (cos(&s)-ry 0 sin(fo , ))(l + r / ) is the complex amplitude transmission 
coefficient of the double mesh out-coupler [12], k = lx/A is the wave number, s is the 
spacing between the meshes, y 0 = (\ + Y 2 + ztan(As))/(l + z(l + I^)tan(As)) is the 

admittance of the double mesh out-coupler, T 7 ={\-y I }l(\+y I ) is the amplitude 
reflectance, and y 2 = y^ + Y^is, the input admittance^ 2], From the previously listed 
equations, it follows that a relationship between Q and 5 can be found, as shown in Figure 
23 (a). The periodic nature of the relationship is the direct result of the periodic nature of 
the interference effects between the meshes, which can be seen in Figure 23 (b) where the 
reflectance, R=\-T, is plotted versus equivalent mesh spacing. 
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Figure 23. (a) Quality factor as a function of mesh spacing. 

(b) Reflectance as a function of mesh spacing. 
Oscillations occur with a period of rik! 2. 
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The two plots in Figure 23 also specify a range of values for the mesh spacing. 
From graphical analysis, 5 is predicted to have a usable range of value between 96pm < s 
< 103pm. Figure 23 (b) shows that if 5 is less than 96pm, the reflectance approaches zero 
and that would not allow net gain in the optical cavity. If 5 is greater than 103 pm, the 
reflectance approaches unity at which point no light is being out-coupled from the cavity. 
Due to the cyclic nature of the interferometer, there are also other values of 5 that yield 
the same optical properties. In this thesis, only the values around 100pm will be 
considered. 

The motivating factor for finding an equivalent mesh spacing for a given Q, flows 
from the prohibitive nature of creating a computational model of the mesh structure. The 
difficulties lie in that the model would need to be on the scale of centimeters with 
openings on the scale of microns. And each opening would need to have enough grid 
points over it to satisfactorily represent the results. This would require large amounts of 
computational time, on the order of weeks per simulation run. In the approximations 
used here, the macroscopic optical properties of the double mesh can be used with a 
degree of confidence in the results. 

C. SIMULATION RESULTS 

With the relationship between Q and 5 established, simulations were conducted 
using the NPS 3D oscillator FEL code and the parameters listed previously, and again 
categorized into two groups based on the electron bunch charge. 

1. lOOpC Bunch Charge Simulations 

Simulations were conducted varying the quality factor over the range 5< Q <100 
in the simulation code, to represent a corresponding mesh spacing, s, as discussed in a 
previous section. The mesh parameters used in the conversion of Q to s were a/g=O. I6 
and g/A= 0.20 which corresponded to a single mesh power reflectivity of 95%. These 
parameters satisfy the condition for long wavelengths. 

Figure 24 shows that as the spacing is increased, the extraction increases to a peak 

and appears to begin to decrease. If the reflectance is too low, there is too little net gain 

in the cavity to exceed the losses. If the reflectance is too high, the fields inside the 
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cavity build up to such a large value that the electrons overbunch, reducing the extraction, 
as discussed in the previous chapter.. The extraction plot also shows that there is a peak 
extraction value of 2.88% when R& 98%which equates to Q=75, and an equivalent 
mesh spacing of s ~ 101/rni. 



Figure 24. The extraction from a double mesh out-coupler with bunch charge 
of 100 pC. The Equivalent Mesh spacing directly corresponds 
to the range of Q : 5<Q<95. 

Figure 25 displays the output diagrams for the peak extraction. It should be noted 
that the mode is close to the fundamental Gaussian as seen in Figure 25 (b). This is 
preferable because the energy produced can then be focused into a smaller spot in an 
application. 
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Figure 25. (a) Simulation output of the evolution of the optical field 

inside the undulator. (b) Optical mode shape 
at the out-coupling mesh surface. 


2. 200 pC Bunch Charge Simulations 

Simulations were also conducted for a bunch charge of 200 pC by varying the 
quality factor over the range of Q values previously stated. Figure 26 illustrates how the 
simulated FEL performed with the increase in bunch charge. The plot shows that the 
peak extraction, 2.9% was approximately the same as the 100 pC case. The optimal 
reflectivity was R= 0.97, which corresponds to a mesh spacing of 5 « 99jum and Q= 35. 
However, the simulated FEL did perform better over a broader range of mesh spacing. 
This can be attributed to the fact that the increase in the bunch charge directly correlates 
to an increase in net gain per pass. 


38 



Figure 26. The extraction from a double mesh out-coupler with bunch charge 
of 200 pC. The green dotted line represents the results 
from the 100 pC case. The Equivalent Mesh spacing directly 
corresponds to the following range of Q : 5<Q<95. 

Figure 23 displays the output diagrams for the peak extraction. It should be noted 
that the mode shape in Figure 27 (b) is once again close to the fundamental Gaussian. 



Figure 27. (a) Simulation output of the evolution of the optical field inside the 

undulator field amplitude in the optical cavity. 

(b) Optical mode shape at the out-coupling mesh surface. 
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3. 


Conclusion 


The double mesh out-coupler simulations showed the FEL had extraction levels 
comparable to hole out-coupling and traditional partially transparent mirrors, as would be 
expected. Unlike hole out-couplers, mesh out-couplers will not disturb the transverse 
mode shape as long as X » g. For larger values of mesh spacing, the 100 pC bunch 
charge may have a larger range of relatively high extraction, whereas the 200 pC bunch 
charge has better performance even at small spacing values. 


40 



VI. FUTURE WORK 


Several THz FEL out-coupling methods were studied in this thesis, including hole 
and double mesh out-coupling. Another method uses zone plate diffractive optical 
elements. Zones plates are surfaces that use diffraction to focus incident light. As seen 
in Figure 28, the zone plate is comprised of concentric rings, or zones, that alternate 
between being transparent and opaque. The transparent zones are spaced to promote 
constructive interference of the light at a specified focal point. To achieve this, the zone 

radii must be determined by r n = JnA / + (n 2 / 1 2 JJ4 , where n= 1, 2, 3, ..., is the integer 

number of the zone, / is the distance to the focus point, and X is the wavelength of light 
being focused [17], 



Figure 28. Zone plate with 5 zones. Blue areas are opaque and 

red areas are transparent. 


Figure 29 shows the diffraction pattern in the plane of a zone plate at the focus 
resulting from five zones. Constructive interference produces a small, central peak at the 
focus, with a faint series of concentric rings due to higher diffraction orders. 
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Figure 29. (a) 2D output intensity of light out-coupled using a zone plate at focal point. 

(b) 3D output intensity of light out-coupled through a zone plate 


This unique method for out-coupling was previously considered to analyze 
whether it was possible to use a zone plate in a resonator [18]. The analysis showed that 
under certain conditions, namely, when the focal length is less than the resonator length, 
a zoned out-coupler would be feasible [18], [19]. 

The NPS 3D FEL code is currently being modified to simulate mirrors with zone 
structures. When this modification is complete, simulations should be conducted in order 
to examine the viability of this method for use in the THz regime. 
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VII. CONCLUSION 


There are many current and potential applications of intense, coherent THz light 
sources, including basic physics experiments, military and homeland security 
applications. The FEL is already a powerful source of coherent light in many parts of the 
electromagnetic spectrum, including THz. From this motivation, this thesis explored 
multiple out-coupling techniques for an FEL operating ~ 100pm. 

Through simulation using the NPS 3D FEL code, hole out-coupling was shown to 
be a viable method for out-coupling in the THz regime. It was seen that the hole had an 
effect on the optical mode in the cavity, and that there was an optimal hole size for 
maximum extraction. This optimal hole size was shown to be affected by changes in the 
bunch charge and the amount of absorption of the mirror material. The hole was also 
shown to affect the transverse mode shape of the beam; when the hole radius was 
increased, alterations in the mode shape increased. 

As suggested from a previous study, beam clipping at the ends of the undulator 
was shown to improve mode quality when using hole out-coupling. An optimal amount 
of clipping was determined for the optimal hole radii studied in this thesis. 

Double mesh out-coupling was also shown to be potentially a good method for 
THz out-coupling. The reflectance of the double mesh out-coupler was shown to 
increase as the spacing between the meshes increased. An optimal reflectance, 
corresponding to an optimal mesh spacing, was determined for different values of bunch 
charge. It was also shown that mesh out-coupling has very little effect on the optical 
mode shape. A double mesh has the additional advantage that it can be adjusted to 
perform well over a broad range of wavelengths, to take advantage of the FEL’s inherent 
tunability. 
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